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Increasing the efficiency of photovoltaic (PV) systems is a pressing issue, 
and several studies have focused on the maximum power point tracking 
(MPPT) techniques to extract the maximum PV output power. Many MPPT 
techniques have been discussed in the last decade, and optimization-based 
MPPT techniques have shown better performance than other MPPT 
techniques. In this study, two optimization techniques, the cuckoo search 
algorithm and particle swarm optimization with changing inertia weight 
techniques, are discussed and applied to a PV system to track the maximum 
power point. The MSX-60 PV module and boost DC-DC converter are used 
in this paper to simulate and model the MPPT system using 
MATLAB/Simulink to show which technique has the best performance 
under various solar irradiation scenarios. In addition, different structures of 
PV arrays such as series-parallel, bridge link, and total cross-tied PV 
structures are simulated to analyze their effect on the efficiency of MPPT 


Photovoltaic systems processes. 
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1. INTRODUCTION 

The burning of fossil fuels has many disastrous effects on our environment as it causes greenhouse 
gas emission, which causes global climate change and the formalization of atmospheric aerosols, increasing 
the global mortality rate [1]. Burning 1 L of fossil fuels can emit approximately 2.9 KG of greenhouse gases 
[2]. Besides, if the use of nonfossil fuel resources is increased by 1%, the emission of carbon dioxide 
decreases by 0.8% [3]. In addition to environmental pollution, fossil fuels are limited, whereas energy 
demands increase as the world population grows; it is predicted to increase by 65% in 2030 [4]. 

Based on these, there is an urgent need to use renewable energy sources, such as hydropower, 
geothermal, solar, biomass, and solar resources. Renewable energy has many advantages: it is free, clean, 
green, and environmentally friendly as it does not cause any greenhouse gas emissions or air pollution [5]. 
Besides the environmental benefits of renewable energies, there are economic benefits, as the use of 
renewable energy can reduce the dependence on imported fossil fuels. In contrast to conventional energy 
resources, renewable energy resources do not deplete. Thus, renewable energy is the best energy resource to 
tackle the environmental and economic challenges faced due to the generation of electricity. 

One of the most favorable resources of renewable energy is solar energy. It can be used directly to 
heat water for households, reducing power consumption by 13% [6]. Electricity can be generated from solar 
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energy using either a direct or an indirect method. The direct method uses photovoltaic (PV) systems, and the 
indirect method uses optical devices that provide steam to propel a turbine and generate electricity [7]. PV 
systems are preferred as they generate no noise, have low preservation expenditure, and can be installed in 
urban areas close to the loads [8], so they are considered as distributed generation (DG) [9]. The most crucial 
advantage of DG is power saving because, during the transmission of power from the power source to the 
consumers, up to 13% of the generated energy can be lost due to the joule effect [10]. 

Despite the advantages of PV systems, the conversion efficiency of PV arrays is extremely low [11] 
and dependent on atmospheric conditions [12], and they have high installation costs [13]. To reduce their cost 
of generating 1 kWh, the PV systems effectiveness must be increased using a dynamic tracking algorithm to 
track the maximum power point (MPP), especially when solar radiation and ambient temperature change. 
There are two types of MPP trackers (MPPTs): mechanical and electronic trackers. Mechanical trackers 
direct the PV panel to track the sun using single or dual-axis trackers [14]. Electronic trackers are preferred 
due to the disadvantages of such tracers, including complexity, low efficiency, and high cost. There are many 
MPPT algorithms that can be categorized into three classes: conventional, soft computing, and optimization- 
based MPPTs. These algorithms differ in convergence speeds, complexity, accuracy, steady-state and 
dynamic efficiencies, number of sensors used in the algorithm, and costs [15]. Many studies have compared 
the listed techniques [16]—[21], and others have improved the existing MPPT algorithms to improve their 
efficiency [22]-[26]. Optimization-based MPPT algorithms have been reported to have a better performance 
than the conventional techniques, especially under partial shading conditions (PSC) [27]-[29], as the 
conventional techniques are a local MPP (LMPP) and cannot track the global MPP (GMPP). 

Many researchers implement particle swarm optimization (PSO) and Cuckoo search algorithm 
(CSA) algorithms to track the MPP, as shown in [30], [31]. It is noted that the PSO algorithm is not very 
efficient in tracking the MPP, so a PSO algorithm with a changing inertia weight strategy is used to enhance 
the algorithm's efficiency and is discussed in this paper. In this paper, PSO with changing inertia weight and 
CSA algorithms are compared to show which one is more efficient in tracking MPP. In addition, different 
structures of PV arrays and their effect on the efficiency of MPPT processes are compared. The rest of this 
paper is arranged as follows: Section 2 describes a PV system; Section 3 presents MPPT techniques; Section 
4 presents the simulation results; Section 5 presents the concluding remarks. The references are listed in the 
final section. 


2. PV SYSTEM 
2.1. PV array modeling 

PV cells are the essential components of a PV array. The most popular single-diode PV cell 
(Figure 1) is composed of a diode in parallel with a photocurrent source in addition to series and shunt 
resistances. These resistances can be neglected for simplicity. 


Figure 1. PV cell model 


The output current of a PV array can be expressed as [32]: 
7 = KpvVopv 
lopy = Npclpg — Nsclrse [exp (Ceereer) -1| (1) 


where Igy, is the PV-array output current, N,- and N,,the number of series and parallel PV cells, 
respectively, Ij, the photogenerated current, and I,,, the reverse saturation current. 


kg, GiGT) (2) 
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Where q, is the electron charge (q, =1.6 xX 10719 C), Kg the Boltzmann's constant (Kg = 1.3805 x 
10-23 J/K), T, the ambient temperature, and I, the ideality factor of the p—n junction. The reverse saturation 
current is calculated using the (3): 


3 
lese = Irs (=) exp|qc-Egn(1/Tra) ona (1/Tq)/1;Ks| (3) 


where T; is the reference temperature, I,,; the reverse saturation current at T;g, E,, the bandgap energy 
(Egp = 1.1 eV). The photogenerated current is calculated as (4): 


Sr 


Ing = toes oP Ag Tra) ) 5 


(4) 
where I,-, is the short-circuit current, K;, the short-circuit current temperature coefficient, s, the solar 
irradiance. The output power is calculated using (5). 


ri KpvVopv 
Pony = lopvVopy = Nocl pg Vopv _ NoclrscVopv (exp (“ere _ 1) (5) 
As MPP occurs at dpDopy/dVopy = 0, the following Equation can be used to calculate the maximum 
output voltage. 


exp (mi) [mn 3] = 6 


Nsc Nsc 


Where v,, is the maximum output voltage. Solving these equations with variation in solar irradiance (s,.), the 
PV-array characteristics can be plotted as shown in Figure 2. Figure 2(a) represents the P-V characteristics, 
and the I-V characteristics are shown in Figure 2(b). 
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Figure 2. PV characteristics under uniform solar irradiance levels: (a) P—V and (b) I-V characteristics 


Figure 2 shows the PV characteristics under uniform solar irradiance, but the PV array may be 
operated under PSC due to buildings, clouds, or dust in natural conditions. Thus, the PV-array characteristics 
change, and several peaks are obtained, as shown in Figure 3. Figure 3(a) represents the P-V characteristics, 
and the I-V characteristics are shown in Figure 3(b). 
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Figure 3. PV characteristics under PSC: (a) P—V and (b) I-V characteristics 
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Because the amount of power generated by a PV cell is small, cells can be connected in parallel or 
series to produce a higher current or voltage, respectively. A PV module is obtained by connecting several 
PV cells. Many modules can be connected to construct a PV array to produce the desired voltage and current. 
PV arrays are connected using three methods [33] (Figure 4). Figure 4(a) represents the series-parallel (SP) 
structure, Figure 4(b) represents the bridge-link (BL) structure, and the total cross-tied (TCT) structure is 
shown in Figure 4(c). 


Figure 4. PV-array structures: (a) SP, (b) BL, and (c) TCT structures 


Modules in the SP structure are connected in series, and each series string is connected in parallel 
with other strings. Each module in the TCT structure is linked in series and parallel with other modules. The 
BL structure is intermediate between the SP and TCT structures. PSC impacts the short circuit current of the 
modules, which affects the output current of the modules at their MPP, resulting in a poor cohesiveness 
between the MPPs of the modules and the MPP of the array. Because the TCT structure has far more parallel 
interlinkages than the SP structure, more current paths are provided to avoid the current decrease due to the 
current decrease in other branches [34]. The TCT structure's parallel interconnections may also decrease the 
probability of turning on bypass diodes, thereby avoiding their relating losses [35]. 


2.2. DC-DC boost converter 

The output voltage of the boost converter is always higher than the input voltage, so it may be called 
a step-up converter. Figure 5 shows a DC-DC boost converter topology. The boost converter consists of a 
voltage source, an inductor, a capacitor, a diode, a switch, and load resistance. 


+ 


I, 
s & Yo 


Figure 5. Boost converter topolog 


The switch (S) can be closed or opened to obtain the desired output voltage. The relationship 
between the input and output voltages is expressed as (7): 


Y= (7) 


where V, is the converter's input voltage, V, is the converter's output voltage, D is the converter's duty cycle. 
The (8) can be used to calculate the average inductor current. 


Vs 


I, = (1-D)2R 


(8) 


Because the output voltage and inductor current depend on the duty cycle, the operational status of the boost 
converter can be controlled by adjusting the duty cycle. An MPPT technique is used to control this duty 
cycle, obtain the desired voltage, and get the boost converter to operate the PV system at the voltage at which 
MPP occurs. 
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3. MPPT TECHNIQUES 

The maximum power of a PV system can be extracted using MPPT techniques irrespective of the 
atmospheric conditions. MPPT uses an electronic system to change the electrical operating point of a PV 
module to deliver maximum power. The schematic diagram of a PV system with the MPPT technique is 
presented in Figure 6. In this section, CSA and PSO algorithms are discussed as they have been proven to be 
more efficient than other techniques 
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Figure 6. Schematic diagram of MPPT 


3.1. Cuckoo search algorithm (CSA) 

Yang and Deb developed CSA in 2009 [36]. It is a metaheuristic algorithm inspired by the breeding 
strategy of cuckoo birds. The CSA algorithm is explained in detail in the flow chart in Figure 7 and the few 
lines below. Cuckoo birds lay their own eggs on other birds' (host birds) nests. To enhance the reproduction 
possibility of their own eggs, cuckoo birds may throw away one of the host birds’ eggs to lay their own egg 
instead of the thrown egg. Cuckoo eggs are often very similar to those of the host birds, except they are slightly 
bigger. Cuckoo eggs are hatched before the host birds' eggs, and once the cuckoo chicks are hatched, they expel 
the host bird’s eggs from the nest to increase their feeding chance. The host bird can discover the cuckoo eggs 
with the possibility p,, where p, € [0,1]. When the host bird discovers cuckoo eggs, it either destroys them or 
relinquishes the nest to build another one. 

Looking for the host nests is similar to the process of searching for food. It is found that Drosophila 
melanogaster and fruit flies search the landscape using a sequence of straightforward flight directions 
interrupted by a rapid 90° shift, which leads to the lévy flight style [37]. In CSA, the eggs in the nest are 
referred to as solutions, and the cuckoo egg is the new solution. The new nests are generated by the lévy 
flight are expressed as (9): 

xD = x + @ @l évy(A), (9) 


L 


where oe is the new solution, aad the current solution, @ the step size (a > 0), and @ the entry wise 


multiplication. | évy(A) can be determined as (10): 
Lévy(A) ~u=l, (10) 


where 1 <A <3. 
The (11) represents a simple scheme of the lévy distribution: 


s=K (<a) (pest — Xi) (11) 


u, v denote the normal distribution curves where: 
u = N(0,07),v = N(0, 0%) (12) 


where o,,, 0, can be defined by (13): 


1/B 
G, = ree Oy = 1 (13) 


Ppa 


where y is the integral gamma function. 
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Initialize the population of (n) host nests and assign the position of 
each nest 


Using Lévy flight get a random nest xi(t) 
Evaluate the fitness of xi(t), F (xi(t}) 


Randomly choose a nest among (n) population, x/(/t) 


Evaluate the fitness of x/(t), F (xj(t)) 


Use x(t) as a new solution 


Replace x//t) by xi(t) using equation (9) 


Minimum error or 


maximum number 
of iterations obtained? 


Find the best solution 


Figure 7. Flowchart of the CSA algorithm 


In MPPT, the converter's duty cycle is represented by the nest position, whereas the fitness 
evaluation function represents the PV-generated power which is represented in (5). 


3.2. Particle swarm optimization (PSO) 

PSO is a swarm-based algorithm introduced by Kennedy and Eberhart [38]. PSO is employed in many 
applications because it is accurate, reliable, easy to implement, requires only a small calculational memory [39]. 
PSO algorithm is explained in detail in the flow chart shown in Figure 8 and the following few lines. PSO was 
inspired by a biological herd of birds seeking food in a specific area. In PSO, the herd is called the “swarm”, 
and the birds are called “particles.” First, the swarm is initialized in the search space with several particles in 
random positions. These particles search for the optimum solution. Each particle's position and velocity are 
updated after each iteration based on (14) and (15): 


si(t +1) = w(t +1) * 5:(t) + City (Doest — Lit) + Cora (Goes — Lilt), (14) 
L(t+1)=1¢)+s5,¢+0), (15) 


where s;(t + 1) is the current velocity; s;(t) the previous velocity; 1;(¢ + 1) is the current position; 1;() is the 
previous position of particle i; ryandr, are random numbers (r,,72 € [0,1]); C, and Cz are the cognitive 
learning rates; w(t + 1) is the inertia weight in the current iteration; Dyes; is the best position achieved by the 
particle; gpe>¢ is the best position achieved by the particles’ neighbors. w can be obtained using (16): 


WMmax—%min 
itermax 


Witer = Vmax — * iter, (16) 
where Wmagx is the maximum inertia weight (Wma, = 0.9), Wmin the minimum inertia weight (Wmin = 
0.4), iter the number of the current iteration, and iter q, the maximum number of iterations. From (13), it is 
evident that the inertia weight is at its maximum value at the beginning of the iterations to accelerate the 
velocity of the particles, which can enhance the global search ability and reach the optimal solution fastly. 
During the iterations, the value of inertia weight is reduced gradually, which results in enhancing the local 
search ability and reaching the optimal solution. After updating each particle's position and velocity, the fitness 
evaluation function is calculated. If the fitness value of the updated particle is better than the memorized 
personal best position (Dyest), Ppest 1S adjusted to make it equal to the position of the current particle, and the 
same steps are iterated for the global best position (Gpes¢). In MPPT, the particle position represents the 
converter's duty cycle, whereas, the fitness evaluation function represents the PV-generated power which is 
represented in (5). 
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Initialize swarm (randomly create every particle with random 
velocity and position) 


Initialize each " pbest’ equal to the current position of each particle, 
Initialize "gbest’ equal to the superior one amongst " pbest'. 


Update iteration counter. 


Update velocity and position for every particle according to equations (14), (15) 
respectively. 


Determine the fitness value for each particle after updating 
If the new fitness value of the particle is superior to the memorized “phest’,set the 
current fitness value as “pbesf" 
If the new fitness value of the particle is superior to the memorized “gbesf’ set the 
current fitness value as “gbest’ 


No 


Minimum error or 
maximum number 
of iterations obtained? 


Yes 


“gbest’ is the optimum solution 


Figure 8. Flowchart of the PSO algorithm 


3.3. Parameter values 

The parameters of metaheuristic algorithms have a significant effect on their performance. As a 
result, the choice of these parameters must be done carefully. The parameter values for both algorithms are 
shown in Table 1. 


Table 1. The parameter values used for both algorithms 


CSA PSO 
Population size 6 Population size 6 
Number of iterations 200 Number of iterations 200 
Pa 0.25 Wmin 0.4 
a 1 Wmax 0.9 
B 1.5 Cy = C2 0.8 
K 0.1 1 =% 1 


4. SIMULATION RESULTS 

The MATLAB/Simulink program was used to build a 3x3 MSX-60 PV array to analyze the 
effectiveness of the MPPT techniques for different PV structures. Four cases were simulated, as shown in 
Figure 9. Case | is represented in Figure 9(a), case 2 is represented in Figure 9(b), case 3 is represented in 
Figure 9(c), and case 4 is represented in Figure 9(d). For every case, PSO and CSA were tested, and the output 
PV power for every case is shown in Table 2. 

Case 1: the standard conditions (STC). In this case, the PV array has been tested with different PV 
structures under STC. The results for this case are represented in Figure 10, Figure 10(a), Figure 10(b), and 
Figure 10(c) represent SP, BL, and TCT structures, respectively. 

Case 2: PSC with bottom shading. In this case, the PV array has been tested under PSC with different 
irradiation levels. The three upper PV modules were irradiated with 1000 W/m7?, the three middle modules 
were irradiated with 900 W/m?and the three bottom modules were irradiated with 700 W/m?. Results for this 
case are represented in Figure 11, Figure 11(a), Figure 11(b), and Figure 11(c) represent SP, BL, and TCT 
structures, respectively. 

Case 3: PSC with side shading. In this case, the PV array has been tested under PSC with different 
irradiation levels. The three top PV modules were irradiated with 900 W/ m?, the three middle modules were 
irradiated with 700 W/m2and the three bottom modules were irradiated with 500 W/m7?. The results for this 
case are represented in Figure 12, Figure 12(a), Figure 12(b), and Figure 12(c) represent SP, BL, and TCT 
structures, respectively. 
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Case 4: In this case, the upper left side of the array was irradiated with 700 W/m2?, and the rest of the 


array was irradiated with 500 W/m?. The results for this case are represented in Figure 13, Figure 13(a), Figure 
13(b), and Figure 13(c) represent SP, BL, and TCT structures, respectively. 


(c) (d) 


Ey frradiation = 1000 w/m? 


(b) 


| irradiation = 900 w/m? 
Ea irradiation = 700 w/m? 


Ey] irradiation = 500 w/m? 


Figure 9. Different solar irradiation scenarios: (a) standard condition (STC), (b) PSC with higher irradiance, 
(c) PSC with lower irradiance, and (d) uniform low solar irradiance 
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Figure 10. Comparison between PSO and CSA performance for Case 1 for different PV structures: (a) SP, (b) 
BL, and (c) TCT 
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Figure 11. Comparison between PSO and CSA performance for case 2 for different PV structures: (a) SP, (b) 
BL, and (c) TCT 
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Figure 12. Comparison between PSO and CSA performance for case 3 for different PV structures: (a) SP, (b) 
BL, and (c) TCT 
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Figure 13. Comparison between PSO and CSA performance for case 4 for different PV structures: (a) SP, (b) 
BL, and (c) TCT 


PSO algorithm has many advantages, such as being accurate, reliable, stable, easy to implement, few 
parameters to be determined, and efficient in global search; also CSA algorithm has many advantages, such 
as requiring only a few parameters to be tuned, and the performance of the CSA algorithm doesn't depend 
heavily on the value of the parameter p,. However, both parameters have serious drawbacks, such as their 
relatively slow convergence speed; also CSA algorithm can easily fall into local optimal solution. 

The simulation results show that the PSO algorithm is more efficient and stable than the CSA 
technique in all tested cases. Obviously, the PSO algorithm has a faster response in most tested cases than the 
CSA algorithm. The simulation results also indicated that the efficiency of the TCT PV structure appears 
especially under SPC. The stability of both algorithms is tested across different runs, each algorithm is 
executed for each tested case, and the variance of solutions is computed using MATLAB's var function, and 
the results are shown in Table 3. From the variance results, it is obvious that both algorithms have good 
stability across different runs, but PSO is more stable than CSA in most of the tested cases. 
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Table 2. The PV output power (W) for Table 3. The variance of output power across 
each case different runs for each case 
Cases Technique SP BL TCT Cases _ Technique SP BL TCT 
Case 1 CSA 519.21 521.50 521.91 Case 1 CSA 0.0037 0.00023 0.002 
PSO 524.16 525.54 528.13 PSO 0.00038 0.00013 0.0018 
Case 2 CSA 365.30 371.47 380.21 Case 2 CSA 4.4532 0.05 0.0430 
PSO 374.58 376.93 386.13 PSO 0.00033 0.00055 0.00042 
Case 3 CSA 254.87 269.9 275.3 Case 3 CSA 0.0063 0.00025 0.0063 
PSO 262.45 278.16 282.69 PSO 0.00018 1.6849 0.00067 
Case 4 CSA 245.67 249.96 261.54 Case 4 CSA 0.0025 0.0039 8.008 
PSO 252.53. 265.28 270.96 PSO 0.00015 0.00037 0.00053 


5. CONCLUSION 

In this study, two optimization-based MPPT techniques were developed, CSA and PSO, and 
simulated them in a PV system using the MATLAB/Simulink program to determine the effectiveness of each 
technique under various atmospheric conditions. The simulation results showed that the PSO technique is 
more efficient and stable than the CSA technique under all tested conditions. The simulation results also 
indicated that varying the PV-array structure can affect the efficiency of the PV system, and the TCT 
structure is more efficient than other PV structures, especially under PSC. However, as explained in the 
following sentences, adapting both algorithms into the PV system presents numerous challenges. Although 
the PSO algorithm tracks the MPPT faster than the CSA algorithm, the PSO algorithm's convergence speed is 
still relatively slow. Although both algorithms achieve excellent stability across different runs in most tested 
cases, they can fall into the local optimal solution, and the result is not assured to be the global optimal 
solution. Finally, both algorithms heavily depend on their parameter values, and changing the parameter 
values will definitely change the result. 
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